Astronomy & Astrophysics manuscript no. final south 


© ESO 2009 


October 6, 2009 





Dense and warm molecular gas in the envelopes and outflows of 

southern low-mass protostars 

T.A. van Kempen^'^, E.F. van Dishoeck''^, M.R. Hogerheijde' , and R. Giisten^ 



o 
o 

(N 

o 
O 



0^ 
O 



o 



'Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands 
2 Center for Astrophysics, 60 Garden Street, MS 78, Cambridge, MA 02138, USA 

Max Planck Institut fiir Extraterrestrische Physik (MPE), Giessenbachstr. 1, 85748 Garching, Germany 
* Max Planck Institut fiir Radioastronomie, Auf dem Htigel 69, D-53121 Bonn, Germany 
e-mail: tvankempenScf a . harvard . edu 



June 2009 



ABSTRACT 



Context. Observations of dense molecular gas lie at the basis of our understanding of the density and temperature structure of proto- 
stellar envelopes and molecular outflows. The Atacama Pathfinder Experiment (APEX) opens up the study of southern (Dec < -35°) 
protostars. 

Aims. We aim to characterize the properties of the protostellar envelope, molecular outflow and surrounding cloud, through observa- 
tions of high excitation molecular lines within a sample of 16 southern sources presumed to be embedded YSOs, including the most 
luminous Class I objects in Corona Australis and Chamaeleon. 

Methods. Observations of submillimeter lines of CO, HCO^ and their isotopologues, both single spectra and small maps (up to 
80" X 80"), were taken with the FLASH and APEX-2a instruments mounted on APEX to trace the gas around the sources. The 
HARP-B instrument on the JCMT was used to map IRAS 15398-3359 in these lines. HCO^ mapping probes the presence of dense 
centrally condensed gas, a characteristic of protostellar envelopes. The rare isotopologues C'*0 and H'^CO^ are also included to de- 
termine the optical depth, column density, and source velocity. The combination of multiple CO transitions, such as 3-2, 4-3 and 7-6, 
allows to constrain outflow properties, in particular the temperature. Archival submillimeter continuum data are used to determine 
envelope masses. 

Results. Eleven of the sixteen sources have associated warm and/or dense (> 10* cm"^^) quiescent gas characteristic of protostellar 
envelopes, or an associated outflow. Using the strength and degree of concentration of the HCO^ 4-3 and CO 4-3 lines as a diagnostic, 
five sources classified as Class I based on their spectral energy distributions are found not to be embedded YSOs. The C'**0 3-2 lines 
show that for none of the sources, foreground cloud layers are present. Strong molecular outflows are found around six sources, with 
outflow forces an order of magnitude higher than for previously studied Class I sources of similar luminosity. 

Conclusions. This study provides a starting point for future ALMA and Herschel surveys by identifying truly embedded southern 
YSOs and determining their larger scale envelope and outflow characteristics. 
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> 1. Introduction 
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During the early stages of low-mass star formation. Young 
Stellar Objects (YSOs) are embedded in cold dark envelopes 
of g as and dust, which absorb the radiation from the central 
star (lLadalll987l lAndreet al]|1993h . This extinction is strong 
enough that low-mass embedded YSOs, or protostars, only emit 
weakly at infrar ed wavelengths (e.g., |j0rgensen et al.l l2005bl 
iGutermuth et al.1 [20081) . Only at later evolutionary phases, 
in which the envelope has been accreted and/or dispersed, 
does emission in the optical and infrar ed (IR) dominate th e 
Spectral Energy Distribution (SED) (e.g. Hartmann et al.ll2005l) . 
Protostars emit the bulk of their radiation at far-IR and sub- 
millimeter wavelengths, both as continuum radiation, produced 
by the cold (T < 30 K) dust, and through molecular line emis- 
sion from the gas-phase species present throughout the proto- 
stellar envelope. Although the bulk of the mass is accreted dur- 
ing the earliest embedded phases, more evolved protostellar en- 
velopes still contain a reservoir of gas and dust that can accrete 
onto the central star and disk system and thus provide the ma- 
terial for disk and planet formation. At the same time, jets and 



winds from the young star interact with the envelope and drive 
molecular outflows which clear the surroundings. Characterizing 
and quantifying all of these different physical components in the 
protostellar stage is still a major observational challenge. 

The protostellar envelopes and molecular outflows can 
be directly observed either through th ermal emission of 
dust at (sub)millimeter wavelengths (e. g. IShirley et alJ l2000l 
iJohnstoneet al J 120011 iNutter et all l2005h or flie line emission 
of molecules. Low frequency molecular emission tra ces the 
cold g as in the protos tellar e nvelopes (e.g., [Ho gerheij de et all 
119981 lJ0rgenseii etal] '2002, Maret et al. 20oJ or molecu- 



lar outflows (e.g. [Snell et al.. ,1990.. Cabrit & Bertoutl 1 19921 



iBachiUer & Tafallalll999l) .' Single-dish observations of dust us 
ing current generation bolometer arrays are able to map large 
areas and image the surroundin gs of pr otostars (e^.^Motte et aLl 
119981 IShirlev etltnimKll IStankeetTanimiellNutter et al.ll2008 
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Through radiative transfer models (e.g. , IShirlev et al.1 l2002l 
|j0rgensen et alj l2002l lYoung et al] l2003l). including info rma- 
tion from shorter wavelengths (e.g.. iHatchell et aT]|20071j ). the 
temperature and density structure of the protostellar enve- 
lope can be constrained, but the continuum data cannot de- 
termine the velocity structure of the infalling envelope, char- 



2 



T.A. van Kempen et al.: Molecular gas associated with southern protostars 



acterize the outflows and their interaction with the surround- 
ings, or disentangle envelope and (foreground) cloud mate- 
rial. Analysis of gas observations in the form of spectra of 
multiple transitions of the same molecule and its isotopo- 
logues provide additional strong constraints on the p hysical 



characteristics of the protostellar envelope (Mangum & Wootten 
1993, Blake et al. 1994, van Dishoeck et al. 1995, Schoier et al 
2002 Maretetal. 2004, J0rgensen et al, 2005b, E vans et al 
2005L Iv an der Tak et al. 2007) and outflowing gas (e.g., 



Cabrit & Bertout 1992, Bontemos et al.'ll996l,'Hogerheiide et alj 
iHatchell et al. 1999, Parise et al. 2006, Hatchefl et al. 

Although many different molecules have been observed in 
protostellar envelopes, only a limited number of species are 
well suited to trace the physical characteristics of all of the 
components of an embedded YSO and its surroundings. For 
example, the use of CH3OH and H2CO is complicated by their 
changing abundances through the envelope, although some 
information can be obtained with c areful analysis and a suffi- 
cient number of observed lines (e.g. Mangum & Wootten 19931 
Blake et al. 1994, van Dishoeck et al. 1995, van der Tak et aQ 
200d iLeurini et alJ l2004l) . The weakness of high-excitation 
CH3OH and H2CO lines in all but a handful of t h e mos t 
luminous Class protostars (e.g. lJ0rgensen et al.l l2005al) 
coupled with a lack of some coUisional rate coefficients 
make these species less suitable tracers for the bulk of the 
low-mass embedded YSOs. In practice, the column density of 
the surrounding cloud (and that of any unrelated foreground 
clouds) is best probed by low excitation optically thin tran- 
sitions of molecules with a low dipole mo ment, such as the 
isotopo logues of CO, e.g., C"*0 2-1 or 3-2 ( van Kempen et al] 
l2009bl) . For molecular outflows the line wings of various 
^^CO transitions have been efficiently used as tracers of 
their properties 



jCabrit & Bertou|_ 
1996', Hogerheiide et al. 19981 IBachiUer & TafaUal 
Hatchell et al. 1999; 



T992I 



Taj). The denser regions of the proto- 
stellar envelope need to be probed with emission lines with a 
high critical density, such as the higher excitation transitions 
from HCO^ with critical densities > 10*' cm"-'. The warm gas 
(T >50 K), which can be present in both the molecular outflow 
and the inner region of the protostellar envelope, can only be 
traced by spectrally res olved high-7 CO transitio ns, which have 
£up > 50 K for 7up > 4 ("Hogerh eiide et al.lll99a) . 

Embedded YSOs are generally identified by their 2-24 mnIR 
slop e with positive values characteristic of Class I objects dLadal 
Il987i) . Over the last several years, it has been found that some of 
these Class I o bjects turn out to be ed^e-on disks or obscured 
sources (e.g. , Luhman & Rieke 1999, Brandner et al.l l2000l 
Pontoppidan et al . 2005, Lahu is et al. 2006, van Ke mpen et alj 
2009bl) . The use of a spectral line map over a small (~2' x 2') 



region around the protostar is an elegant solutio n to disentan- 
gle th e contributions of the different components jBoogert et alj 
l2002j) . In particular, the spatial distribution of the HCO+ 4-3 line 
has proven to be an excellent diagnostic of such 'false' embed- 
ded sources: in the case of Ophiuchus, about 60% of the sources 
with C lass I or flat SE D slopes turned out not to be embedded 
YSOs jvan Kempen et al. 2009b). 

With the development of the Atacama Large 
Millimeter/Submillimeter Array (ALMA) on Chajnantor, 
Chile, surveys of sources in southern star-forming clouds 
will be undertaken at high resolution (< 1") at a wide range 
of frequencies, and are expected to reveal much about the 
inner structure of protostars and their circumstellar disks. The 
Herschel Space Observatory will also target a large number 



of low-mass YSOs across the sky. Both sets of observations 
rely heavily on complementary large aperture ground-based 
single-dish observations for planning and interpretation. Apart 
from providing the necessary information about the structure on 
larger scales, single-dish studies of low-mass star formation in 
the southern sky will also put the results obtained from studies 
on the northern hemisphere in perspective. 

Due to a lack of sub-millimeter telescopes at high dry 
sites in the southern hemisphere, high frequency data on south- 
ern YSOs are still limited. The Swedish ESO Submillimeter 
Telescope (SEST) operated up to the 345 GHz window, but 
for only limited periods of the year As a result, southern 
clouds, such as Chamaeleon, Corona Autralis, Vela or Lupus, 
are much less studied than their counterparts in the north- 
ern sky, such as Taurus and Perseus, where such observa- 
tions have been readily availa ble for ov er a decade (e.g 



Hogerheiide et al. 1997, Motte et all Il998l [Hogerheiid e et al 



1998, Johnst one et al.u200a ,j0rgensen et al.ll200Z .Nutter et al 
12005; 2008). The Atacama Pathfinder Experiment (APEX) 
dGiisten et alj |2006) Q has opened up access to the atmospheric 
windows in the 200-1400 jum wavelength regime over the entire 
southern sky. 

Of the southern star-forming regions not or only poorly vis- 
ible with northern telescopes two regions have proven espe- 
cially interesting. The Chamaeleon I cloud (D=130 pc, Dec = 
-77°), observed with IRAS and the Infrared Space Observatory 
(ISO) and included in the Spitzer Space Telescope guaranteed 
time and 'Cores to disks' (c2d) Legacy programs jPersi et al 



^200di lEvans et all I2003L iDamianovet 311120071 iLuhman et al 
l2008h contains some embedded sourc e s, especially around 
the Cederblad region ('Persi et al.' '200d', ' Lehtinen et all 1200 iL 
jielloche et al. 2006, Hiramatsu et al. 200^ Corona AustraUs 
is a ne arby star-forming region (D=170 pc, see iKnude & H0gl 
(Il998l) ; Dec - -36°), well -known for the central Coronet clus- 
ter near the R CrA star (|Loren| [UtI iTavlor & StorevI 11984 
IWilking et al]|1986l lBrownlll987h . Large-scale C'^O 1-0 maps 
show that it contains about 50 M© of gas and dust (Harju et aTl 
I99 3I). Many surveys have been undertaken in t he IR (e.g., 
Wil king et alJ [T986t [19971 lOlofsson et al] [1991 iNisini et al.1 
2005^ but only a few studies mapped this region at s ubmillime- 
ter wavelengths dChini et al.ll2003LlNutter et alj|2005h . Although 
Corona Australis has only a few protostars with rising infrared 
SEDs, the cloud does contains some of the most luminous low- 
mass protostars in the neighborhood of the Sun {D < 200 pc) 
with luminosities of up to 20 L©. 

In this paper we present observations of submillimeter lines 
of CO and HCO^ of a sample of 16 embedded sources in the 
southern sky, with a focus on the Chamaeleon I and Corona 
Australis clouds, to identify basic parameters such as column 
density, presence and influence of outflowing material, presence 
of warm and dense gas and the influence of the immediate sur- 
roundings, in preparation for Herschel and ALMA surveys or 
in-depth high resolution interferometric observations. In section 
§ 2 we present the observations, for which the results are given 
in § 3 with a distinction between the clouds and isolated sources. 
Both the single spectra (§ 3.1) and maps (§ 3.2) are presented. 
In § 4 we perform the analysis of the observations, making use 



' This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX) with programs E-77.C-0217, E- 
77.C-4010 and E-78.C-0576. APEX is a collaboration between the 
Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 
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Source 


RA 


Dec 


Dist. 


Lum. 


Ref." 


Notes" 




(h m s) 
[J2000] 


(d m s) 
[J2000] 


(pc) 








Chamaeleon I 


Ced 110 IRS 4 


11 06 47.0 


-77 22 32.4 


150" 


1 


1,2,3 


ISO-Cha 84 


Ced 110 IRS 6 


1 1 07 09.6 


-77 23 04.3 


150 


0.6 


1,3 


ISO-Cha 92 


Cha IRS 6a 


1 1 07 09.9 


-77 23 06.3 


150 


0.8 


4 


ISO-Cha 92+ 


ChaIRN 


11 08 37.1 


-77 43 51.0 


150 


5 


3 


ISO-Cha 150 


Cha INa 2 


11 09 36.6 


-76 33 39.0 


150 


0.6 


5 




Corona Australis 


RCrA IRS 5 


19 01 48.0 


-36 57 21.6 


no- 


2 


6,14 


SMM5 


HH 100* 


19 01 50.7 


-36 58 10 


no 


15 


6,14 


SMM3 


RCrA IRS 7A 


19 01 55.2 


-36 57 21.0 


170 


3.3'' 


7,14 


SMA 2, SMM IC 


RCrA IRS 7B 


19 01 56.2 


-36 57 27.0 


170 


17'' 


7,14 


SMA 1, SMM IB 


RCrA TS 3.5 


19 02 07.1 


-36 53 26.2 


170 








CrA IRAS 32 


19 02 58.7 


-37 07 34.5 


170 


3.4 


14 


ISO-CrA 182, SMM 8 


Isolated 


HH46 


08 25 43.8 


-51 00 35.6 


450-' 


16 


8,9,10 


BHR36 


IRAS 07178-4429 


07 19 21.8 


-44 34 55.1 






10,11 


CGl/BHR 17 


IRAS 12496-7650 


12 53 17.2 


-77 07 10.6 


200« 


50 


12 


DK Cha, Cham II 


IRAS 13546-3941 


13 57 42.2 


-39 56 21.0 


550'' 




10,11 


CG12/BHR 92 


IRAS 15398-3359 


15 43 01.3 


-34 09 15.0 


130' 


0.92 


2,13 


B228, Lupus 



' ' Refe rences: IjLuhma n et al.| (|2008D; 2. i FroebricH j2005l) ; 3. iPersi et alj (I200C ) ; 4. iPersi et al.l (120011): 5. iPersi et alj (Il999h: 6. Nisini etal.1 
('2005'); 7. Groppi et al. (2007); 8. Norieaa-Crespo et al. (2004); 9. Velus amv et al.l ( |2007|) : lO. lBourke et all jl995h : ll. lSantos et al.l ( ll998l) : 12. 
[van Kempen et al. (2006); 13. E vans et a l. ( 2005J; 14. Nutter et al. (2005.) 

Observations of HH 100 were erronously pointed at coordinates: RA 19h01m49.1s, Dec -36d58ml6.0s. See text for further discussion. 

Other names and cloud location. 

'' The luminosities of the two RCrA IRS 7 sources are not well determined, due to confusion of the total integrated flux at far-IR wavelengths. The 

t otal luminosity of both sources (20 Lq) was used in the characterisation. 

'■^ude&H0g(1998) 

f Heathcote et al. (1996) 

^ Hughes & Hartigan ( 199"1) 

''iMaheswar et al. (20041) 

' iMurphv etal..a986i) 



of archival submillimeter continuum data, with the final conclu- 
sions given in § 5. 



2. Technical information 

2.1. Source list 

Sources were selected from a sample of sources with rising in- 
frared SEDs (Class I) observed with the InfraRed Spectrograph 

i IRS) on Spitzer in the scope of the c2d legacy program 
Evans et aTl 12003!) and with the Very Large Telescope (VLT) 
of the European Southern Observatory (Pontoppidan et al. 2003h 
(Table 1). At the time of selection of the Spitzer and VLT sources 
in 2000, the list included the bulk of the known southern low- 
mass Class I YSOs with mid-infrared fluxes of at least 100 mJy. 
Most of these sources are located in the Chamaeleon and Corona 
Australis clouds, supplemented by a few other isolated Class I 
sources found initially by IRAS. Within the Chamael eon I cloud, 
most sources are located in the Cederblad region dPersi et al] 
[2000, Hiramatsu et al. 2007). Cha IRS 6a is located 6" away 
from Ced 110 IRS 6 (-hal f an APEX bea m at 460 GHz). Ced 
110 IRS 4 was studied by IBelloche et al.l (12006) which detect 
an outflow with an axis of <30°, which makes this source a 
prime target in Chamaeleon. The embedded YSOs in the Corona 
Australis region are mos tly located in the small regi on around R 
CrA, called the Coronet dWilking et al . 1986; 1 99llNutter et alj 
I2005h . but two sources are located in the R CrA B region. RCrA 



IRS 7A and B are believed to be two embedded objects in a 
wide binary configuration with a separation of 30" dNutter et al.l 
I2005L ISch5ier et aLllIOO^ . These are deeply embedded as little 
to no IR counterparts were found (Groppi et al. 2007). HH 100 
is in a region that is dominated by the RCrA A complex, 
but is an embedded source in itself according to iNutter et"al] 
(120051) . CrA IRAS 32 is located in R CrA B region, but has 
little to none known surrounding material. HH 46 is a Class I 
protostar famous for its well-char acteriz ed outflow, pointed to 
flie south-we st (Chernin & Mass onlll 99 iL [Heathcote et alJI 19961 
Stanke et aliri999> .Noriega-Crespo et al.'7004, Velusamy et alj 
2007h . It is extensively discussed in|yan Kempen et al. (2009a|. 
IRAS 12496-7650, also known as DK Cha, is'one of the bright- 
est protostars in t he southern hemisphere a nd a potential driving 
source of HH 54 dvan Kempen et alj|2006h. ISO-LWS observa- 
tions were detected with ISO (*Gianni ni et al ]|199 9'). but not con- 
firmed by Ivan Kempen et al . (2006). IRAS 07178-4429, IRAS 
13546-3941 and IRAS 15398-3359, afl very interesting embed- 
ded sources i i i the s o uthern hemisphere, a re also included (e.g . , 
iBourke et all [1991 IShirlev etalJ [lOOOl iHaikala et all llooa) . 
Table [T| lists the sources of our sample and their properties. 
Liiminosities are adopte d from the c2d lists dEva ns et al. 20031 ; 
l2009l) and the VLT lists dPontoppidan et al.ll2003b . 

During the reduction, it was discovered that the original po- 
sition of HH 100 was incorrect and differed by ~20" from the 
IR position. Throughout this paper, the incorrect position will 
be refen-ed to as 'HH 100-off', while 'HH 100' will refer to the 
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actual source. Single-pixel spectra for HH 100-ofF are shown in 
some of the figures, as they are useful as probes of the cloud ma- 
terial and outflows in the RCrA region, whereas the true HH 100 
source position is included in the maps. 



Table 2. Overview of the APEX observations 



Line 




Instrument 


Cloud" 


Note 




3-2 


APEX-2a 


Cha, CrA, Iso. 






map 


APEX-2a 


Various 


80" X 80" 




4-3 


FLASH-I 


Cha, Iso 






map 


FLASH-I 


Cha, Iso 


40" X 40" 




7-6 


FLASH-II 


Cha, Iso 






map 


FLASH-II 


Cha, Iso 


40" X 40" 




3-2 


APEX-2a 


Cha, CrA, Iso. 




HCO+ 


4-3 


APEX-2a 


Cha, CrA, Iso. 






map 


APEX-2a 


Various 


80" X 80" 


H'3CO+ 


4-3 


APEX-2a 


Cha, CrA, Iso. 





" Cha = Chamaeleon, CrA = Corona Australis and Iso. = Isolated sam- 
ple 



2.2. Observations 

Observations of the s ources in T able [T] were carried using the 
APEX-2a (345 GHz, iRisacher et al. 2006) and the First Light 
APEX Submillim eter Heterodyne (460/800 GHz, denoted as 
FLASH hereafter. iHevminck et al.ll2006) ) instruments mounted 
on APEX between August 2005 and July 2006. FLASH allows 
for the simultaneous observation of molecular lines in the 460 
GHz and 800 GHz atmospheric windows. 



IRAS 15398-3359 






{ IRAS 


13546-3941 


IRAS 0717B-4429 




HH 46 f 








-20 -10 10 20 

Velocity (km/s) 



-20 -10 10 20 

Velocity (km/s) 



Fig. 3. I^CO 7=3-2 (left) and C'^O 3-2 (right) spectra for the 
isolated sample. 

Using APEX-2a, observations were taken of the '^CO 3-2 
(345.796 GHz) , C'^O 3-2 (329.331 GHz), HCO+ 4-3 (356.734 
GHz) and U^^CO^ 4-3 (346.998 GHz) molecular Unes of the 
entire sample at the source position with the exception of 
Cha IRS 6a in '^CO. The '^CO 4-3 (461.041 GHz) and 7-6 



IHAS ISSde-SO-ll 

10 ^^^^A^^^f^^^Aj^^r 



IRAS 12496-7650 




IRAS 13546-13941 



IRAS 12496-7650 



IRAS 07178-4429 




-10 10 

Velocity (km/s) 



-10 10 
Velocity (km/s) 



Fig. 4. HCO+ 7=4-3 (left) and H^^CO^ 4-3 (right) spectra for 
the isolated sample. 




Velocity (km/s) 



Velocity (km/s) 



Velocity (km/s) 



Fig. 5. '^CO 7=4—3 (left) and 7-6 (middle) spectra of the sources 
in Chamaeleon, as well as the spectra of HH 46 and IRAS 
15398-3359 taken with FLASH (right). 



(806.652 GHz) transitions were observed with FLASH for the 
Chamaeleon I sample as well as some of the isolated sample 
(Ced 110 IRS 4, Ced 110 IRS 6, Cha INa 2, Cha IRN, HH46, 
IRAS 12496-7650 and IRAS 15398-3359). The single specti-a of 
CO 3-2, C' ^O 3-2, CO 4-3 and 7-6 for IRAS 12496-7650 can 
be found in V an Kempen et aTl (l2006h . but the integrated intensi- 
ties are included here for co mpleteness. Simila rly, HH 46 is dis- 
cussed more extensively in van Kempen et al.l ('2009a). Typical 
beam sizes of APEX ai-e 18", 14" and 8" at 345, 460 and 805 
GHz with respective main beam efliciencies of 0.73, 0.6 and 0.43 
(iGiisten et al.ll2006l) . 

For both instruments, new Fast Fourier Transform 
Spectrometer (FFTS) units wer e used as back- ends, with 
over 16,000 channels available jKlein et al.l l2006h . allowing 
flexible observations up to a resolution of 60 kHz (0.05 km 
s ' at 345 GHz). Pointing was checked using line pointing on 
nearby point sources and found to be within 3" for the APEX-2a 
instrument. The pointing of FLASH was not as accurate with 
excursions up to 6" for these observations. However, it is not 
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-10 10 
Velocity (km/s) 




-10 10 
Velocity (km/s) 



-20 -10 10 
Velocity (km/s) 



10 10 

Velocity (km/s) 



Fig. 1. ^^CO 3-2 spectra of the sources in Chamaeleon (far left) and Corona Australis (middle left) and C'^O 3-2 spectra of the 
sources in Chamaeleon (middle right) and Corona AustraHs (far right). Spectra in this and subsequent figures are offset for clarity. 




Ced 110 IRS 6 



j^^-^^^mV^^^^ 




Ced 110 IRS 4 



Cha INa 2 
Cha IRN 



Ced 110 IRS ( 



RCrA IRS 7A 




RCrA TS 3,5 



HH 100 - off 



CrA IRAS 32 



-20 -10 10 20 -20 -10 10 20 -20 -10 10 20 -20 -10 10 20 



Velocity (km/s) 



Velocity (km/s) 



Velocity (km/s) 



Velocity (km/s) 



Fig. 2. HCO^ 4-3 spectra of the sources in Chamaeleon (far left) and Corona Australis (middle left), as well as H'^CO^ 4-3 spectra 
of the sources in Chamaeleon (middle right) and Corona Australis (far right). 



expected that this has significant influences on the results, as 
most spectra were extracted from the peak intensities of the 
maps. Calibration was done on hot and cold loads as well as 
on the sky. Total calibration errors are of the order of 20%. 
Reference positions of 1,800" to 5,000" in azimuth were 
selected. For RCrA IRS 7 and 5, emission in the CO 3-2 line 
at this position was subsequently found resulting in artificial 
absorption in the spectra, but this did not exceed 5 % of the 
source emission and is less than the calibration uncertainty. 

Observations at APEX were taken under normal to good 
weather conditions with the precipitable water vapor ranging 
from 0.4 to 1 .0 mm. Typical system temperatures and integration 
times were 600 and 1000 K and 5 and 4 minutes for APEX-2a 
and FLASH-I, respectively. All spectra at a given position were 



summed and a Unear baseline was subtracted within the CLASS 
package. The APEX-2a observations were smoothed to a spec- 
tral resolution of 0.2 km s ' . The RMS was found to be 0. 15-0.2 
K in 0.2 km s ' channels for HCO"" 4-3, '^CO 3-2 and C"^0 

3- 2. H'^CO^ spectra were re-binned to 0.5 km s ' for deriving 
the (upper limits on) integrated intensity. The RMS for the CO 

4- 3 observations is 0.35 K and for CO 7-6 1.4 K in channels of 
0.3 km s ', significantly higher than the APEX-2a observations 
because of the higher system temperatures. 

In addition to single spectra, small (up to 80" x 80") fully 
sampled maps were taken using the APEX-2a instrument (see 
Tabled of '^CO 3-2 and HCO+ 4-3 for HH 46, Ced 110 IRS 4, 
CrA IRAS 32, RCrA IRS 7 and HH 100 . With FLASH, similar 
maps were taken of the Cha I sample, HH 46 and IRAS 12496- 
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7650 in '^CO 4-3 and 7-6, but not of the RCrA sample. All 
maps were taken in on the fly (OTF) mode and subsequently 
rebinned to 6" . The RMS in the maps is a factor of 5 or more 
higher than that in the single spectra. For example, for the CO 
4-3 maps taken with FLASH the noise levels are ~0.9 K in a 0.7 
km s"' channel, compared with 0.35 K in a 0.3 km s"' for the 
single position spectra. 

For IRAS 15398-3359, maps were also taken with the 
16-pixel heterodyne array receiver HARP-B mounted on the 
James Clerk Maxwell Telescope (JCMTfl in the CO 3-2, C"*0 
3-2 and HCO^ 4-3 lines. The high spectral resolution mode 
of 0.05 km s"' available with the ACSIS back-end was used 
to disentangle foreground material from outflowing material. 
Spectra were subsequently binned to 0. 15 km s ' . The HARP-B 
pixels have typical single side-band system temperatures of 
300-350 K. The 16 receivers are arranged in a 4x4 pattern, 
separated by 30". This gives a total footprint of 2' with a spatial 
resolution of 15", the beam of the JCMT at 345 GHz. The 
2' X 2' fields were mapped using the specifically designed jiggle 
mode HARP-^^- Calibration uncertainty is estimated at about 
20% and is dominated by absolute flux uncertainties. Pointing 
was checked every two hours and was generally found to be 
within 2-3". The map was re-sampled with a pixel size of 
5", which is significantly larger than the pointing error. The 
main-beam efficiency was taken to be 0.67. The JCMT data for 
IRAS 15398-3359 are presented in § 3.3 

All APEX data were reduced using the CLASS reduction 
package 0. For the JCMT data the STARLINK package GAIA 
and the CLASS package were used. 



3. Results 

3.1. APEX Single spectra at source position 

Spectra are shown in Figs. [T]-|4] (all spectra taken with APEX- 
2a), and Fig. |5] (spectra taken with FLASH). The integrated in- 
tensities and peak temperatures at the source positions are given 
in Tables [3] (CO and C'**0 observations with APEX-2a and 
FLASH observations) and|4](HCO+ and H'^CO^ with APEX- 
2a). Source velocities are determined with an accuracy of 0.1 
km s~' by fitting gaussian profiles to the optically thin lines of 
H'^CO"^ (preferred if detected) or C'^O. The results are pre- 
sented in Table |4] Table |3] also includes the width of the C"*0 
3-2 line, if detected. These are generally narrow, of order 1 km 
s-i. 

In general the 3-2 and 4-3 lines of '^CO cannot be fitted 
with single gaussians, because of self-absorption at or around 
the source velocity Vlsr- This absorption is especially deep in 
RCrA IRS 7 A and 7B with t >> 1. Due to this absorption, the 
total integrated line strength of '^CO spectra is not useful in the 
analysis of the quiescent molecular gas. Only for a few sources, 
e.g. IRAS 07178-4429, can the '-CO 3-2 be fitted with a single 
gaussian. Interestingly, the self-absorption is absent in all three 
detected CO 7-6 lines in Ced 110 IRS 4, Cha INa 2 and Cha 

^ The James Clerk Maxwell Telescope is operated by The Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. Data 
were obtained with programs M06BN1 1, M07BN09 and M08AN05 

^ See JCMT website http://www.jcmt.jach.hawaii.edu/ 
CLASS is part o f the GIL DAS reduction package. See 
[http://www.iram.fr/JRAMFTt/GILDA S for more information. 



IRN at the S/N level of the observations. All C"*0 3-2 lines 
can be yfitted with single gaussians, indicating that there are no 
unrelated (foreground) clouds along the line of sights. The C'^O 
3-2 lines are much stronger in the R CrA cloud than in the other 
clouds. 

The majority of the sources show strong line wings, a sign 
of a bipolar outflow. Especially the sources in RCrA have promi- 
nent emission up to 15 km s ' from the line center. 

The HCO^ lines, which can be used to trace the dense gas, 
cover the widest range of antenna temperatures and include a 
variety of line profiles. As two extreme cases, no HCO^ 4-3 is 
seen toward Cha INa2 down to 0.2 K, whereas the spectrum for 
RCrA IRS 7 A peaks at rMB = 18 K. Note that the brightest spec- 
tra in Corona Australis all show self-absorption, indicating that 
the HCO^ 4-3 is optically thick at the positions of the envelope. 
Ced 110 IRS 4, RCrA IRS 7 A and 7B and CrA IRAS 32 have 
line profiles suggesting dense outflowing gas. The H'^CO^ 4-3 
line is only detected for 5 sources, although the detection of Ced 
110 IRS 4 is only 3. 4o-. 

3.2. APEX maps 

The maps of the total integrated intensities of HCO^ 4-3 can be 
found in Fig. |6]and|2] alongside outflow maps of '^CO 3-2. In 
the outflow maps, the integrated spectrum is split between red- 
shifted (dashed contour), quiescent (not shown) and blue-shifted 
(solid contour). The three components are separated by ± 1 .5 km 
s"' from the systemic velocity (see Table 4). Different methods 
for separating the outflow and quiescent components were tested 
but these did not change the results within the overall uncertain- 
ties. The integrated intensity maps of CO 4-3 are shown in Fig. 
[8] These are comparable to the CO 3-2 maps, but not with the 
HCO+ 4-3 maps. 

The HCO+ 4-3 and CO 3-2 maps for Ced 110 IRS 4 and 
CrA IRAS 32 (see Fig. |6l) show spatially resolved cores, cen- 
tered on the known IR position. For both sources, outflow emis- 
sion is seen originating at the source position. HH 100 (Fig. 
|6ll shows no HCO^ 4-3 core, but does have outflowing gas in 
the CO 3-2 map, originating at the IR position (marked with 
a X). There is strong HCO^ emission north of HHIOO, which 
does not coiTespond to any source in the iNutter et al I (l2005b 
SCUBA maps. At this moment, this discrepancy cannot be ex- 
plained. The maps around RCrA IRS 7A and 7B (Fig.|7] sources 
marked with X and centered on RCrA IRS 7A) show a single 
elongated core in the HCO^ map that is up to 60" along the ma- 
jor axis. Both sources lie firmly within the core. The '^CO 3-2 
outflow is centered on RCrA IRS 7A, suggesting it is the driving 
source. HH 46 (Fig.|7j shows a nearly circular core in HCO^, 
but a strongly outflow-dominated emission profile in CO 3-2, 
with much stronger re d-shifted than blue - shifted emission. For 
more information, see Ivan Kempen et alj (l2009al) . The CO 4-3 
integrated intensity maps show elongated spatial profiles for Ced 
1 10 IRS 4, Ced 1 10 IRS 6 and Cha IRN, but no peak for Cha INa 
2. IRAS 12496-7650 shows a round core, dominated by outflow 
emission (see § 4.3.3) 

3.3. JCMT map of IRAS 15398-3359 

Figure|9]shows the spectrally integrated intensity maps of CO 3- 
2, HCO+ 4-3 and C"*0 3-2 around the position of IRAS 15398- 
3359. The maps clearly show a single core around the central 
position, which is completely spatially resolved. There is little 
to no C'**0 or HCO^ emission outside of the core, while the 
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Fig. 6. HCO+ 7=4-3 (left) and CO 3-2 (right) maps of Ced 1 10 IRS 4 (top), CrA IRAS 32 (middle) and HH 100 (bottom). The maps 
are 80" x 80" in size. HCO^ 4-3 is integrated over the entke line. CO 3-2 shows the red- (dashed lines) and blue-shifted (solid 
lines) parts of the emission, derived by integrating the line wings of each spectral line (velocity difference greater than 1.5 km s~' 
w.r.t. the systemic velocity). Contour levels for the CO 3-2 map are at 3cr,6cr, 9cr, ... with cr equal to the noise levels of 0.3 K km 
s~'. The HCO^ contours are plotted at 10%, 20%, 30%, ... of the maximum intensity given in Table|4] The 10% level corresponds 
to la- for Ced 1 10 IRS4, 4o- for CrA IRAS 32 and 9cr for HH 100. 

The APEX beam is shown in the lower right image. The IR position of HH 100 is marked with a 'X', while the pointed position HHlOO-off is 
marked with a plus sign. 



CO 3-2 has emission lines at around 10 K in the surrounding 
cloud. The C"^0 map, which lacks outflow emission, is some- 
what shifted to the south from the other two maps. The spectra 
at the central source position taken from these maps can be found 
in figure [10] The CO spectra are dominated by the blue-shifted 



outflow. The integrated intensity of HCO^ 4-3 can be found in 
Table|4] The independently calibrated CO 3-2 spectra taken with 
JCMT HARP-B differ from the spectrum taken with APEX-2a 
by only 3% in integrated intensity, well within the stated calibra- 
tion uncertainties of both telescopes. 
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Fig. 7. HCO+ 7=4-3 {left) and CO 3-2 (right) maps of RCrA IRS 7 {top) and HH 46 {bottom). The maps are 80" x 80" in size. 
HCO^ 4-3 is integrated over the entire Hne. CO 3-2 shows the red- {dashed lines) and blue-shifted {solid lines) emission, derived 
by integrating the Hne wings (velocity difference greater than 1.5 km s"' w.r.t. the systemic velocity) of each spectral line. Contour 
levels for the CO 3-2 map are at 3cr,6o-, 9cr, ... with cr equal to 0.3 K km s The HCO+ contours are plotted at 10%, 20%, 30%, 
... of the maximum intensity given in Table |4] The 10% level corresponds to 25cr for RCrA IRS 7, and 5cr for HH 46. The APEX 
beam is shown in the lower right image. Both RCrA IRS 7A and IRS 7B are shown with X in the CO map, but the map is centered 
on RCrA IRS 7A. 



4. Analysis 

4.1. Envelope properties 

Assuming the (sub-)mm emission is dominated by the cold dust 
in the protostellar envelop e, its mass can be d e rived from dust 
continuum observations (I Shirlev et all l2000l ^J0rgensen et al 



2002 ). Using the SCUBA Legacy archive (IDi Francesc o et al 



20081) . en velope masses were derived for YSOs in Corona 
AustraUs. iHenning et al] (Il993l) observed fluxes at 1.3 mm for 
a large sample in Chamaeleon using the SEST, which are used 
here to constrain the masses of the very southern sources. All 
masses were derived using the relation 



M,o,^SyD^/KyBy{Ti) 



(1) 



assuming an isothermal sphere of 20 K, a gas-to-dust ratio of 1 00 
and a dust emissivity Ky at 850 fim of 0.02 cm ^ g"'(gas-i-dust) 
from lOssenkopf & HenningI (column 5. Il994l) . To scale to 1.3 
mm, the dust emissivity was assumed t o scale with frequency 
as Ky PC v'^ following the results from 'Ossenkop f & HenningI 
(11994 . The mass of the HH 46 envelope is estimated at ~5 
Mq using new LABOCA observatio ns and radiative transfer 
modelHng dvan Kempen et all l2009al) . For IRAS 12496-7650, 



LABOCA data of Cha II was used (Nefs, van Kempen & van 
Dishoeck, in prep) The radii Rcore over which Sy is ca l culate d 
are taken to be the core radii from ' Di Francesco et alJ (l2008h . 
Such radii are not the FWHM of the cores, but a re quite close to 
the r adius where the temperature reaches 10 K (|j0rgensen et al.l 
12002) . No mass could be derived for Cha IRS 6a, RCrA TS 
3.5, IRAS 07178-4429, Cha INa2 and IRAS 13546-3941 as they 
were not detected or observed. The estimated envelope masses 
are given in Table|5] They range from 0.04 to 6.3 M©. For Corona 
Australis, our mas ses are a factor of three lower than those by 
iNutter et alJ (l2005l) . Differences are due to the difference in as- 
sumed dust temperature of 12 vs 20 K, as well as a different dust 
emissivity Ky at 850 yum. 

With the assumption that the H^^CO^ 4-3 and C'^O 3-2 
emission is optically thin, one can estimate the column densi- 
ties of these molecules in the protostellar envelope. In this cal- 
culation, an excitation temperature of 20 K is assumed, based 
on the average temperature within the protostellar envelope, for 
both the envelope and cloud material. Table |5] gives the column 
densities of H2, C"^0 and H'^CO^, where the H2 column den- 
sity is derived from the C'**0 column de nsity assuming a H2/CO 
ratio of lO-* and a '<'0/"*0 ratio of 550 dWilson & Roodl[T994b . 
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Fig. 8. CO 7=4-3 integrated intensity maps taken with FLASH of Ced 1 10 IRS 4, Ced 110 IRS 6 and Cha INa2 (top row), and Cha 
IRN, HH 46 and IRAS 12496-7650 (bottom row). Note that the noise levels in these maps are higher than in the spectra given in 
Fig S Source positions are marked with 'X'. The APEX beam is shown in the lower right image. The lines are drawn to guide the 
eye. 

Table 3. '^CO and C"*0 results'" 



Source 


'^CO 3-2 




C"*0 3-2 






CO 4-3 




CO 7-6 






/ TuBdV 




/ TuBdV 




AV 


/ TMBdV 


T'mb 




T'mb 




(Kkm s-') 


(K) 


(Kkms-') 


(K) 


(km s-') 


(Kkms-') 


(K) 


(Kkms-') 


(K) 


Chamaeleon 


Ced 110 IRS 4 


30.7 


7.5 


3.7 


2.1 


1.5 


29.8 


7.5 


24.8 


6.9 


Ced 110 IRS 6 


20.0 


9.3 


1.6 


1.4 


1.1 


9.7 


4.9 


8.8 


2.5 


Cha IRS 6a 






1.2 


1.1 


1.0 










Cha IRN 


11.1 


3.5 


2.3 


1.8 


1.2 


6.5 


3.0 


14.0 


5.2 


Cha INa 2 


11.3 


6.7 


1.3 


1.6 


0.8 


10.0 


4.7 


<1.5 




Corona Australis 


CrA IRAS 32 


51.6 


II.O 


4.6'- 


4.4 


1.0 










HH 100-off 


102.6 


24.7 


6.3 


5.2 


1.1 










RCrA IRS 5 


144.2 


29.9 


11.4 


8.7 


1.2 










RCrA IRS 7A 


406.0 


46.4 


22.5 


8.9 


2.2 










RCrA IRS 7B 


332.7 


40.4 


23.0 


10.0 


2.2 










RCrA TS 3.5 


16.8 


6.4 


<0.3 














Isolated 


HH 46'' 


82.5 


194 


3.2 


3.3 


0.9 


70.5 


14.9 


46.5 


8.6 


IRAS 07178-4429 


24.0 


114 


1.7 


2.3 


0.7 










IRAS 12496-7650* 


92.8 


25 


9.2 


1.9 


1.1 


90.0 


23 


43.7 


19.2 


IRAS 13546-3941 


16.3 


11.4 


1.8 


2.3 


0.7 










IRAS 15398-3359 


25.8 


9.0 


0.5 


1.2 


0.4 


19.9 


5.6 


45.3 


8.5 



" Integrated intensities are across the entire line profile. 
* From van Kempen et al. (2006|). 
^ CO 7-6 from van Kempen et al. ( 2009a). 
Calibration uncertainties estimated at 20% dominate the u ncertainty of the integ rated intensity. 

This is a significant discrepancy from the value reported in lSchoier et alj(l2006h. which determ ined this intensity to be 6.5 K km s-'. A possible 
explanation can be the inaccurate calibration during science verification of lSchoier et al.l ( l2006h . 



The H2 column densities range from a few times 10^^ cm-^ in 
Chamaeleon to a few times 10^^ cm"^ in Corona Australis, all 
in a 18" beam. However, freeze-out of CO onto the interstel- 
lar grains make the adopted ratio of H2 over CO a lower limit. In 
very cold regions (T <15 K), abundances can be as much as two 
orders of magnitude lower than the assumed CO/H2 abundance 



of 10 . For sources where continuu m data at 850 /im are avail- 
able from Di Fran cesco et al.l (1200 8). the column density was in- 
dependently calculated from the dust, assuming a dust temper- 
ature of 20 K. These numbers are in general a factor of 2 to 
3 higher consistent with some freeze-out of CO, but for some 
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Table 4. HCO+ and H^CO^ results. 



Source 


HCO+ 4- 


3 


H"CO+ 


4-3 


Vlsr" 






Tub 


/ ^MBdV* 


Tub 






(Kkms-') 


(K) 


(Kkms-') 


(K) 


(km s-') 


Chamaeleon 


Ced 110 IRS 4 


4.0 


3.1 


0.38 


0.25 


4.3 


Ced 110 IRS 6 


0.5 


0.4 


<0.06 




4.8 


Cha IRS 6a 


<0.2 




<0.06 




4.6 


Cha IRN 


0.87 


0.95 


<0.06 




4.5 


Cha INa 2 


<0.2 




<0.06 




4.8 


Corona Australis 


CrA IRAS 32 


4.1 


2.9 


0.72 


0.53 


5.6 


HH 100-off 


4.8 


3.2 


<0.06 




5.9 


RCrA IRS 5 


14.2 


7.0 


0.96 


0.86 


5.7 


RCrA IRS 7A 


63.7 


17.4 


3.8 


1.6 


5.5 


RCrA IRS 7B 


49.0 


18.7 


3.8 


2.2 


5.7 


RCrA TS 3.5 


<0.2 




<0.06 




5.9 


Isolated 


HH46 


8.4 


5.4 


1.1 


1.0 


5.3 


IRAS 07178-4429 


1.0 


0.6 


<0.08 




3.4 


IRAS 12496-7650 


1.9 


1.5 


<0.06 




1.9 


IRAS 13546-3941 


<0.12 




<0.06 




-5.7 


IRAS 15398-3359'- 


5.9 


4.6 






5.2 



" The rest velocity Vlsr was determined by fitting gaussians to either the H'^'CO^ 4—3 (preferred) or C"*0 3-2 lines (see Table[3}- 
* Calibration uncertainties estimated at 20% dominate the uncertainty of the integrated intensity. 
^ Data from HARP-B map. 

Table 5. Overview of the envelope properties". 



Source 






lope A'(H2)" 


iV(C'»0) 


A'(H''CO+) 


TC03-2 


ThCO+4-3 




(") 


(Mo) 


(10^2 cm-2) 


(10'5 cm-2) 


(10'^ cm--) 






Chamaeleon 


Ced 110 IRS 4 


20 


0.07 


5.0/- 


8.3 


1.4 


181 


6.7 


Ced 110 IRS 6 


20 


0.04 


3.3/- 


5.5 




90 


<13 


Cha IRS 6a 






0.26/- 


0.43 








Cha IRN 


20 


0.17 


4.3/- 


7.2 




400 


<5.2 


Cha INa 2 






3.8/- 


6.3 




150 




Corona Australis" 


CrA IRAS 32 


44 


0.8 


10.5/21.2* 


17.5 


3.0 


257 


16 


HH 100-off 


45 


1.6 


12.4/35.2* 


20.7 




130 


<0.9 


RCrA IRS 5 


43 


1.7 


20.7/24.5* 


35.1 


4.8 


189 


10 


RCrA IRS 7A 


48 


6.3" 


21.2/30.5* 


35.3 


8.9 


117 


7.7 


RCrA IRS 7B 


48 


6.3" 


23.8/30.5* 


39.7 


12.3 


156 


10 


RCrA TS 3.5 






-/- 










Isolated 


HH 46" 


60 


5.1 


8.1/18.2* 


13.5 


6.1 


103 


16 


IRAS 07178-4429 






5.5/- 


9.2 




123.9 


<8 


IRAS 12496-7650 


25 


0.4 


3.0/- 


5.0 




25 


<3.1 


IRAS 13546-3941 






5.5/- 


9.2 




123 




IRAS 15398-3359 


25 


0.5 


1.9/- 


3.2 




79 





" The column density of H2 is derived from that of C'**0 assuming H2/C"*0=5.5xlO*. If 850 fim data are available in lDi Francesco et all d2008h . 
column densities based on the dust emission are also computed. These numbers are indicated as the second entry with an asterisk. Both column 
densities are within the same 18 " beam. 

* Mass estimates in the RCrA region may be severely overestimated due to presence of extended cloud material. 

" RCrA IRS 7A and 7B cannot be distinguished in the sub-millimeter continuum maps; the value refers to the core containing both sources. 
'' For a thorough discussion on the mass of HH 46, see .van Kempen et al. ( .2009a) . 



sources the column densities derived from the dust approach 
those obtained from C'^O. 

The isotopologue observations of C"*0 and H^^CO^ also al- 
low to estimate the average optical depth r in the CO 3-2 and 
HCO^ 4-3 lines (Table |5]l following the formula of 

/(CO)//(C'*0) = (1 - e-^)/(l - e-^'^^")) (2) 



with the isotopologue ratios for both C'^O and H'-'CO^ taken 
from I Wilson & RoodI ([T994.) . For I, the peak temperatures from 
Tables 3 and 4 were taken to avoid contributions from outflow 
material. The optical depths range from 25 to a few hundred 
in the CO 3-2 and from about 0.9 to 16 for the HCO+ 4-3. 
HCO^ 4-3 is still quite optically thick, even though it is asso- 
ciated with just the envelope. One error in the estimates from 
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Fig. 9. The spectrally integrated intensities of CO 3-2 left, HCO+ 4-3 middle and C'**0 3-2 Wg/if observed toward IRAS 15398- 
3359 with HARP-B on the JCMT. Spectra were integrated between velocities of -5 and 15 km s The source itself is at 5.2 km 
s ' . The cores are all slightly resolved, compared to the 15" beam of the JCMT. 
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Fig. 10. The spectra of CO 3-2 top, HCO+ 4-3 middle and C'^O 
3-2 bottom observed with HARP-B on the JCMT at the central 
source position. 



the optical depths, may be that the C'^O 3-2 and H'^CO+ them- 
selves are not optically thin, but have optical depths of a few. In 
that case, the optical depth of the main isotope would lower as 
well. However, observations of even rarer isotopologues, such as 
C'^O or HC"*0+ would be needed. 

Recent C"^0 3-2 data of the Ophiuchus cloud show that at 
many locations throughout L 1688, multiple foreground layers 
are present providing ad ditional reddening of embedded sources 
dvan Kempen et al.l200 9b). In our sample, all C'**0 3-2 lines can 
be fitted with single Gaussians indicating that no fore-ground 
layers are present in any of the observed clouds, in contrast with 
Ophiuchus. 



4.2. Embedded or not? 

Analysis of the physical structure of several sources classified as 
Class I based on their IR spectral slope has shown that some of 
them are in fact not embedded sources, but rathe r edge-on disks, 
such as CRBR 2422.8 (Brandneret al. 2000, P ontoppidan et all 
120051 ), or IRS 46 ( Lahuis et al..,2006,) . Others turn out to be ob- 
scured (backgroun d) sources or occasiona lly even disks in front 
of a dense core (van Kempen et al. 2009b). Theoretical mod- 
els from Robitaille et al. (2006) and Crapsi et al. (2008) indeed 
show that edge-on disks can masquerade as embedded sources 
in their infrared SED due to the increased reddening. 

Following iRobitaille etall ( |2006h and ICrapsi et all ( l2008l) . 
the border between a truly embedded so-called 'Stage 1' source 
and a 'Stage 2' star + disk system is put at accretion rates of 
10"^ Mq yr"', which in the context of these models corresponds 
to an envelope mass of ~0. 1 M© or a H2 column of (2 - 6) x 10^' 
cm"^. Because of potential confusion with surrounding (fore- 
ground) clouds and disks, one cannot simply use the masses de- 
rived from the continuum data but one must look in more de- 
tail at the properties of a source. One of the key additional cri- 
teria for identification as a truly embedded source is the con- 
centration of warm dense gas around the source. The concentra- 
tion is determined by the comparison of the peak intensity in a 
map with the amount of extended emission. Such a parameter 
has been extensively discussed and successfully applied in the 
analysis of submilUmeter continuum observations from SCUBA 
(iJo hnstone et al.ll200lL IWalawender et al.ll2005l |j0rgensen et al.l 
2007^ For example, they have shown that pre-stellar cores have 
much more extended and less concentrated emission than proto- 
stellar envelopes. 

A recent survey of Class I sources in the Ophiuchus L 1688 
region has demonstrated that such a method can be expanded 
to molecular line mapping and be used as a powerful tool for 
identifying truly embedded protostars from edge-on disks and 
obscured sources (Ivan Kempen et aLll2009bh . HCO+ 4-3, a high 
density tracer, proved to be particularly useful. The recipe for 
proper identification consists of three ingredients. First the peak 
intensity of HCO^ line should be brighter than 0.4 K, a facto r 
of 4 more than is commonly seen for disks (iThi et al.ll2004l) . 
Second, the concentration parameter of HCO^ 4-3 should be 
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Table 6. HCO^ 4-3 and CO 4-3 core radii and concentrations 



Source 


Rhco* (") 


Chco+ 


^CO(4-3) (") 


CcO(4-3) 


Emb.* 


Chamaeleon 


Ced 110 IRS 4 


21 


0.77 


20 


0.69 


y 


Ced 110 IRS 6 






16 


0.74 


y 


Cha IRS 6a 










n" 


ChalRN 






12 


0.72 


y 


Cha INa 2 








n.c." 


n'' 


Corona Australis 


CrA IRAS 32 


15 


0.78 






y 


HH 100 




n.c." 






y^ 


RCrA IRS 5 










r 


RCrA IRS 7A 


25 


0.68 






y 


RCrA IRS 7B 


25 


0.66 






y 


RCrA TS 3.5 










n'' 


Isolated 


HH46 


15 


0.79 






y 


IRAS 07178-4429 




n.c.° 






n 


IRAS 12496-7650 












IRAS 13546-3941 




n.c." 






n 


IRAS 15398-3359 


22 


0.9 






y 



" n.c. = no core associated with source 

' Our assessment whether the source is embedded or not 

Based on presence of outflow, see §4.2, 4.3 and Fig|6] 
'' no HCO^ detected down to the limit. 

Strong HCO+ seen. 



higher than 0.6. Third, the HCO^ should preferably be extended 
within the single dish beam. For non-resolved HCO^ emission, 
other information, such as whether or not the source is extended 
in the submillimeter continuum or whether it has a (compact) 
outflow, should be taken into account. 

The concentration of HCO^ is defined as: 

_ l.UB'Suco^ 

LHCO+ - 1 Tj— = , (Jj 

with 5hco+ the HCO^ emission integrated spatially (K km s ') 
and spectrally over the entire envelope with radius Robs (in arc- 
seconds), and Tint - J TMsdV the HCO^ emission integrated 
spectrally over the central beam B of 18" alone (in K km s 
The observed core radius is taken to be the FWHM of the spa- 
tial distribution of the integrated HCO^ intensity (in arcseconds). 
These values of Robs, tabulated in Table|6] are often comparable 
to the core radii found from the continuum data. Although less 
reliable due to a lower critical density and potential contributions 
from outflows, the quiescent part of the CO 4—3 emission, which 
is associated with the central regions of protostellar envelopes, 
can serve as a diagnostic of the concentration of the warm dense 
gas in the inner envelope as well. Note that any outflow wings 
have to be removed. 

Table |6] lists the derived concentration parameters for our 
sources. Of the 12 sources for which such data exist, 6 have con- 
centrations of ~0.7 and higher, as well as strong HCO^ fluxes 
above 0.5 K km s"' associated with the sources, characteristic 
of embedded sources. Three sources (HH 100, IRAS 07178- 
4429 and IRAS 13546-3941) have no associated HCO+ core 
and thus no concentration. Similarly, Cha INa2 shows no signs 
of centrally concentrated CO 4-3 emission. Cha IRS 6a, Cha 
INa2, RCrA TS 3.5, IRAS 07178-4429 and IRAS 13546-3941 
all lack HCO^ 4 -3 emission, even below the limit adopted from 
iThi et al.l(l2004b . and are thus very likely not embedded YSOs. 



HH 100 is an exception. Although no concentration could 
be calculated for this source as it is not peaking on the HH 100 
source position, there is strong extended HCO^ emission with at 
least 7cr, and in some areas up to 50cr. A clear outflow is seen 
in the CO emission in Fig. |6l thus indicating that HH 100 is an 
embedded source. This is confirmed from results from Spitzer- 
IRAC and MIPS (Peterson, priv. comm.). However, the lack of 
concentration from an envelope is puzzling due to the low optical 
depth derived for this line at HHIOO-OFF. 

The RCrA IRS 7 binary sys t em wa s classified as two em- 
bedded sources by iGroppi et al.l (l2007h . The elongated HCO^ 
spatial emission profile indeed suggests that two envelopes are 
located close to each other The concentrations of 0.68 and 0.66 
(close to 0.7) together with the strong outflow confirm that both 
sources are embedded. Similar arguments are used to classify 
RCrA IRS5 and IRAS 12496-7650 as embedded YSOs. In sum- 
mary, using our line diagnostics, we have confirmed that eleven 
of our sixteen Class I sources are truly embedded YSOs, and five 
sources very likely not. 

4.3. Outflows 
4.3.1. Temperatures 

If we assume that these line wings fill the beam and are trac- 
ing the dominant part of the outflowing gas mass (also known as 
the swept up gas), then the ratios of emission in the line wings 
of the CO 3-2, 4-3 and 7-6 profiles at the central position in 
Table [T] can be compared to model line ratios from Appendix 
C of Ivan der Tak et alj (l2007i) . derived using statistical equilib- 
rium calculations within an escape probability formalism. The 
line wings are assumed to be optically thin and filling factors 
to be the same within the different beams. A quick look at the 
line wings of CO 3-2 and C'^O 3-2 in the sources in Corona 
Australis provides an upper limit to the optical depth of the 
CO line wing of about 3.1 close to the quiescent part of the 
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Table 7. Temperature constraints of CO line wings ratios''. 



Source 


CO 3-2/4-3 


CO 4-3/7-6 




,„(K) 


7'outflow (K) 












Scenario 1" 


Scenario 2* 




Red"^ 


Blue'^ 


Red 


Blue 


Red 


Blue 


Red Blue 


Ced 110 IRS 4 


1.0 


1.2 


3.5 




50 


40 


170 160 


Cedll01RS6 


2.2 


1.7 






<30 


40 


50 60 


HH46 


1.0 


1.6-2.0 


2.0 




50 


<40 


180 80 


IRAS 12496 




0.8 




>4 




50 


<200 


IRAS 15398 




0.7 




0.8 




100 


>200 



" Lines thermally excited (n > «„) 

* Lines subthermally excited (« assumed to be lO'' cm"^) 

Red and Blue outflow contributions are calculated by excluding the central 3 km s"' emission. 

Errors on the temperatures are 25%, except for temperatures above 150 K, where it is estimated to be at 40%. 



Table 8. Outflow parameters from the CO 3-2 mapping. 



Source 


i 


Mass" 




P b 






^CO 








(Mo) 


(km s-') 


(AU) 


(yr) 


(Mo 


(Mo yr-' 


(Lo) 














yr-') 


km s-') 




Red Lobe 


Ced 110 IRS 4 


15 


1.4e-2 


4.0 


2.4e3 


2.8e3 


4.2e-6 


4.6e-5 


7e-5 


CrA IRAS 32 


45 


6.6e-2 


6.3 


9.3e3 


7.1e3 


9.3e-6 


1.3e-4 


5.5e-4 


HH 100 


60 


3.6e-l 


14.7 


>6.1e3 


>2.0e3 


<1.5e-4 


<4.9e-3 


<7.2e-2 


RCrA IRS 7 


45 


l.leO 


23.0 


>le4 


>2.1e3 


<4.3e-4 


<2.2e-2 


<5.0e-l 


HH46 


35 


4.3e-l 


9.5 


>2.5e4 


>1.2e4 


<3.2e-5 


<4.4e-4 


<3.1e-3 


IRAS 15398-3359 


75 


2.9e-2 


8.2 


2.9e3 


4.6e3 


6.4e-7 


6.4e-6 


1.6e-5 


Blue Lobe 


Ced 110 IRS 4 


15 


5.8e-3 


-4.2 












CrA IRAS 32 


45 


1.7e-2 


-5.4 


4.6e3 


4.0e3 


4.3e-6 


5e-5 


1.8e-4 


HH 100 


60 


4.7e-2 


-9.7 


4.6e3 


2.2e3 


2.1e-5 


4.4e-4 


2.9e-4 


RCrA IRS 7 


45 


1.2e0 


-15.8 


>5e3 


>1.5e3 


<7.9e-4 


<2.7e-2 


<2.9e-l 


HH46 


35 


1.4e-l 


-7.4 


1.4e4 


9.0e3 


1.5e-5 


1.7e-4 


9.3e-4 


IRAS 15398-3359 


75 


9.2e-2 


0.7 


3.9e3 


4.1e4 


2.2e-6 


3.3e-5 


1.2e-4 



" Corrected for inclination using the average correction factors of lCabrit & Bertou3 il990h . 
* Not corrected for inclination. 

Dynamical time scale : = i?/Vmax 

Mass outflow rate : M = M/ti. 
' Outflow force : Fco = MVl,JR 
f Kinetic luminosity : Lki„ = Q.SMVl^^^^jR 



line. At more extreme velocities, the opacity is probably much 
lower and indeed optically thin. Optically thin line wings were 
also shown to be a good assumption for the case of HH 46 by 
Ivan Kempen et alj (l2009ah . Kinetic temperatures can then be de- 
rived if outflow emission is detected in two of these three lines. 
Note that an observed ratio leads to a degenerate solution of 
the kinetic temperature and local density and not a unique so- 
lution of the temperature. Thus, some estimate of the local den- 
sity must be available. Densities in the outer envelope can be 
estimated from radiative transfer modelling of the dust contin- 
uum emission and SED dlvezic & Eh tzur 1997). H2 densities of 
a few times 10^ cm-^ are derived by'j0rgensen et al. (2002) for 
distances of 1500-3000 AU (coiTesponding to 10-20" at the dis- 
tances of Chamaeleon and Corona Australis), but these drop to 
a few times 10'* cm"^ at a few thousand AU. These densities are 
assumed to be comparable to the densities in the swept-up gas at 
these distances. 

At densities above the critical density, ncr=5xlO"^ cm-^ for 
CO 3-2, the gas is thermally excited and a temperature can be 
inferred directly. For lower densities, the gas is subthermally ex- 



cited and higher kinetic temperatures are needed to produce sim- 
ilar ratios of the line wings. To estimate the temperatures, two 
scenarios are presented. In the first, it is assumed that the gas 
is thermally excited (n > n^) since densities in the envelope are 
near to, or higher than, the CO 3-2 critical density. This is appli- 
cable for all sources with the exception of HH 46, for which the 
central beam extends to about 4,500 AU in radius. The kinetic 
temperatures determined in this way are lower limits. A second 
scenario is considered by assuming that the densities are of order 
10"^ cm"^, closer to typical densities of the surrounding cloud or 
envelopes at larger radii. In the first scenario, outflow temper- 
atures are on the order of 50 K. An exception is IRAS 15398- 
3359, which appears to be much warmer (100 K) from both the 
CO 3-2/4-3 ratio as well as the 4-3/7-6 ratio in the blue out- 
flow. Since we lack CO 4-3 and/or 7-6 data of the red outflow, 
no temperature could be determined for the red outflow. If den- 
sities are lower (scenario 2), derived kinetic temperatures are on 
the order of 150 K, although for HH 46 (blue) and Ced 1 10 IRS 6 
(both red and blue) the outflow temperatures are consistent with 
~80 K. Errors on the outflow temperatures are generally on the 
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order of 25 to 40 %, due to the sensitive dependency on density, 
as well as errors on the estimate of the CO line ratios. 

Bow-shock driven shell models by Hatchell et al.l d 19991) pre- 
dict temperatures of the swept-up gas by calculating the energy 
balance between the heating due to the kinetic energy of the ex- 
panding momentum-conserving shell and the line cooling of the 
gas. For parameters typical of the Class L 483 source (L=9 Lq, 
Menv= 4.4 Mq), outflow tcmpcraturcs of 50 to 100 K are found 
in the inner few thousand AU from the star. Along the bulk of 
the outflow axis, temperatures are modelled to be 100-150 K. 
Higher temperatures in excess of a few hundred to a thousand 
K are predicted only near the bow shocks. The models depend 
on the local conditions such as dust cooling, density distribu- 
tion and jet velocity, but most parameters have an influence of 
the order of a few tens of K on the temperature along most of 
the outflow. The result s in Table [T] agree w ell with the tempera- 
ture predictions from Hatchell et al.l (11999 ) and suggest that jet 
velocities >100 km s ' are not present, with the exception of 
IRAS 15398-3359. 

Although we assumed the line wi ngs to b e optically thin, 
this is not nec essarily the case (e.g., ICabrit & Bertout 1991 
iHogerheiideet al.iri998l) . A more in-depth discussion of the ef- 
fects of optical depth and density on the kinetic temperatures of 
the swept -up gas in the molecular ou tflow of HH 46 IRS is pre- 
sented in Ivan Kempen et akl (l2009al) . A lower line ratio for the 
same temperature can also be caused by a higher optical depth 
in one of the line wings. Similarly a higher density will lead 
to lower inferred tem peratures for the same observed line ratio. 
iHatchellet 311 ( 119991) also show that there may be a temperature 
gradient away from the connecting line between bow shock and 
origin at the protostar However, this gradient is often diluted 
out and may be observable within a single beam. The different 
temperatures of two different ratios may be responsible for this. 

4.3.2. Other outflow parameters 

For the six sources where outflow emission has been mapped 
in '^CO 3-2, it is possible to determine the spatial outflow 
parameters, such as outflow mass, maximum velocity and ex- 
tent. Subsequently the dynamical time scales and average mass 
outflow rate, the outflow force (or momentum flux) and ki- 
netic luminosity can be derived following a recipe detailed in 
iHogerheiide et al.l (Il998h with the assumption that the emis- 
sion at the more extreme velocities is optically thin and that 
the material has kinetic temperatures of 50-80 K. For IRAS 
15398-3359 a temperature of 100 K is adopted. The inclination 
of the outflow with the plane of the sky influences the results 
and the derived values need to be corrected fo r their inclination 
using an average of the model predictions in ICabrit & BertoutI 
( 19901). For HH 46 the inclination has been con strained to 35° 
feeipurth & Heaflicoteiri991l iMicono et aT]|199 8'). For the other 
sources, an inclination is estimated from the distance between 
the major and minor axes of the maximum contour levels in Figs. 
|6]and|7](see Tabled. 

The results in Table [8] show that most outflows have masses 
of 0.01-0.1 Mo, with average mass loss rates between 2x10"^ 
to 4.4x10""* Mq yr"', and outflow forces Fqo of 6.4x10"*' to 
2.7x10"^ Mo yr"' km s"'. The strongest outflow is RCrA IRS 
7, which produces the highest outflow force and the most kinetic 
luminosity. This is mainly caused by the fact that emission is 
detected out to AVmax as high as 23 km s"', compared with 4 
to 9 km s"' for the other sources, almost an order of magnitude 
difference. Note that for HH 100 (red only), HH 46 (red only) 
and RCrA IRS 7 (red and blue), the quoted outflow radius is a 



lower limit, due to insufficient coverage. The resulting outflow 
force and kinetic luminosity are thus upper limits. For HH 100 
and and RCrA IRS 7 subtracting Gaussians with widths of 3 km 
s"' changes the integrated intensity by 20%. All other param- 
eters depend linearly on the mass M that is derived from the 
integrated intensity. 

IHogerheiide et alj ([1998) surveyed the outflows of Class I 
sources in the Taurus cloud. There, most sources produce weaker 
outflows than found in Table [Sj with outflow forces of the order 
of a few times 10"^ M© yr"' km s"', with a few sources having 
values as high a few times 10""* M© yr"' km s"' . A recent survey 
of the outflows present in the Perseus cloud shows Fco values 
of a fe w times 10"^ Mp yr"' k m s"' to a few times 10"^ Mq yr"' 
km s~' (iHatcheU et al.l20()7al) . However, these numbers were not 
corrected for inclination, which can account for at least a factor 
of 5 increase. The observed outflow forces in our sample are 
all among the upper range o f these values, especia lly HH 100 
and RCrA IRS 7. Figure 5 of iHatchell eTall (l2007al) shows that 
most Class I outflows have lower outflow forces than the Class 
sources. 

Most flo ws observed from the no rthern hemisphere were first 
discussed in lBontemps et alj (Il996h using CO 2-1 observations 
with relatively large beams. Outflow forces range from a few 
times 10"^ to lO""* Mg y r"' km s~\ corre cted for inclination 
and optical depth. In lBontemps et al.l ( 11996 ). a relation between 
the outflow forces, bolometric luminosities and envelope masses 
was empirically derived, log(Fco) - -5.6-1- 0.91og(Lboi) and 
log(Fco) = -4.15 + 1.1 log(Me„v). 

Figure (TTj shows a comparison bet ween the observed out - 
flows of Table |8] and the relations from lBontemps et al.l ( Il996l) . 
Although the outflows in our sample are all associated with Class 
I sources, they are exceptionally strong outflows for their lumi- 
nosity, but not for their envelope mass. Outflow forces are one 
to two orders of magnitude higher than expected from the above 
equation with luminosity. This suggests that the outflow strength 
as measured through the swept-up material says more about the 
surroundings (i.e., the amount of matter that can be swept up) 
than about the intrinsic source and outflow properties. 

4.3.3. The outflow of IRAS 12496-7650 

An outflow wa s identified from the line profiles in 
Ivan Kempen et al.l (l2006l) and i t was concluded th at the 
emission as seen with ISO-LWS dGiannini et al. I I2OOTI) has to 
originate on scales larger than 8" and not from the inner regions 
of a protostellar envelope around IRAS 12496-7650. Figure [T2] 
reveals that the detected blue-shifted outflow is located only 
in the central 15" and is peaking strongly on the source. The 
lack of red-shifted emission reinforces the conclusion that the 
outflowing gas is viewed almost perfectly face-on and that our 
line of sight is straight down the outflow cone. 

4.3.4. The outflow of IRAS 15398-3359 

The line profiles of IRAS 15398-3359 have a clear blue outflow 
component, with red-shifted emission notseen until further off 
source. If only the line wings are mapped. Figure [13] shows that 
the outflow of IRAS 15398-3359 is smafl and extends only -25" 
across the sky. As discussed above, the blue part of the outflow 
is exceptionally warm, with significant outflow emission seen in 
the CO 7-6 line. Although it is smaller in extent, the outflow 
does not stand out w.r.t to the others in terms of outflow proper- 
ties as measured from their CO 3-2 maps. 
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Fig. 11. CO outflow force Fco versus the luminosity, Lboi (left), and versus the envelope mass, M^w (middle). The ri^ht figure 
show s Fco/L\,oi versus Menv/i^^p which should be free of most luminosity and distance effects following Fig. 7 of iBontemps et alJ 
([1996). Red outflows are shown with diamonds and blue outflows with triangles. The relations between Fco, iboi and Menv found 
by Bontemps et al. are shown with dashed lines. 



IRAS 12496-7650 CO 4-3 




10 -10 

A RA ["] 



20 




Fig, 12. CO 7=4-3 map of IRAS 12496-7650. The blue and red- 
shifted outflow emission are shown in solid and dashed contours 
respectively in steps of 10% of the maximum of 34.8 K km s"'. 

5. Conclusions 

We present observations of CO, HCO^ and their isotopologues, 
ranging in excitation from CO 3-2 to CO 7-6 of a sample of 
southern embedded sources to probe the molecular gas con- 
tent in both the protostellar envelopes and molecular outflows in 
preparation for future ALMA and Herschel surveys. The main 
conclusions are the following: 

- HCO^ 4-3 and CO 4-3 integrated intensities and concen- 
trations, combined with information on the presence of out- 
flows, confirm the presence of warm dense quiescent gas as- 
sociated with 11 of our 16 sources. 

- RCrA TS 3.5, Cha IRS 6a, Cha INa2, IRAS 07178-4429 and 
IRAS 13546-3941 are likely not embedded YSOs due to the 
lack of HCO^ 4-3 emission and/or the lack of central con- 
centi-ation in HCO+ 4-3. 



Fig. 13. CO 7=3-2 map of IRAS 15398-3359. The blue and red- 
shifted outflow emission are shown in blue and red contours re- 
spectively in steps of 10% of the maximum of 8.1 K km s-1. 
The grey-scale and dashed contour lines in the background are 
the 850 fim continuum. 



- The swept-up outflow gas has temperatures of the order of 
50-100 K, as measured from the ratios of the CO line profile 
wings, with the values depending on the adopted ambient 
densities. These values are comparable to the temperatures 
predicted by the heating model of Hatchell et al. ( 1999), with 
the exception of IRAS 15398-3359, which may be unusually 
warm. 

- The outflows of 6 of our truly embedded sources — Ced 110 
IRS 4, CrA IRAS 32, RCrA IRS 7 A, HH 100, HH 46,IRAS 
15398-3359 — were characterized using CO spectral maps. 
The outflows all have exceptionally strong outflow forces, 
almost two orders of magnitude higher than expected from 
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their l uminosities following the relation of iBontemps et al.l 
([T996I) . 

- Neither Chamaeleon nor Corona Australis have foreground 
layers as found in Ophiuchus L 1688. All C'**0 3-2 spectra 
can be fitted with single gaussians. 

Future observations using ALMA and Herschel will be able 
to probe the molecular emission associated with these YSOs at 
higher spatial resolution and higher frequencies. Comparison 
with single dish data, both continuum surveys and spectral line 
mapping, will be essential in analyzing the results obtained with 
future facilities. 
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